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Z.  INTRODUCTION 


The  Plasma  Research  Group  of  the  University  of  California 
at  Davis  continued  the  theoretical  work  of 
microwave-atmospheric  interactions,  in  close  coordination  with 
personnel  at  the  Naval  Research  Laboratory  and  elsewhere. 
This  report  covers  the  period  September  1,  1981  to  August  31, 
1982.  We  have  made  advances  in  two  areas:  (1)  expanding  the 
previous  survey  into  new  regimes;  and  (2)  simulations  of  re¬ 
cent  NRL  focused  microwave  experiments.  ^^The  main  results  of 
our  invesitgation  are  described  in  the  following.  Typical  hy¬ 
drodynamic  calculations  show  that  an  ionization  front  is  ra- 
pldly  formed  which  moves  toward  the  microwave  source  and  con¬ 
sequently  decouples  the  microwaves  from  the  original  ioniza¬ 
tion  region.  By  focusing  the  microwaves  or  using  a  reflector, 
ionization  can  be  confined  to  localized  regions  where  the  mi¬ 
crowave  strength  is  high  enough  to  cause  breakdown  even  though 
the  incoming  microwaves  are  below  threshold.  In  a  strongly 
collisional  atmosphere,  it  is  found  that  the  cutoff  plasma 
density,  which  roughly  equals  the  collisionless  cutoff  density 
multiplied  by  the  collisional ity,  is  much  higher  than  the  cal¬ 
culated  maximum  density.  This  results  in  high  absorption  that 
increases  with  microwave  power,  decreases  with  atmospheric 
pressure,  and  is  quite  independent  of  other  parameters.  In  a 
weakly  collisional  atmosphere,  breakdown  easily  creates  a 
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plasma  density  higher  than  the  cutoff  density  and  causes  re¬ 
flection.  It  is  found  that  the  reflection  decreases  with  the 
collisionality  of  the  system  and  is  quite  independent  of  the 
microwave  strength.  In  general,  the  microwave  field  strength 
In  the  Ionization  regions  is  attenuated  to  the  breakdown  value 
at  steady  state,  and  the  resulting  electron  temperature  Is 
about  2  eV,  independent  of  the  incident  microwave  flux. 
Limitations  of  the  calculations  are  due  to  the  availability  of 
experimental  data  for  the  rate  coefficients,  but  comparison  to 
results  from  recent  focused  microwave  experiments  shows  excel¬ 
lent  agreement. 

In  the  past,  most  studies  and  measurements  were 
done 1  near  the  breakdown  thresholds.  Recently  experiments 
have  been  performed  ,  for  microwave  power  higher  than  the 
threshold,  to  study  the  behavior  of  the  plasma  after  break¬ 
down.  Our  investigation  will  cover  the  overall  scope  of  the 
phenomena  and  will  be  described  in  more  detail  in  the  follow¬ 
ing. 


If  we  neglect  free  electron  diffusion,  experimental 
data3  show  that  the  breakdown  threshold  is  roughly  given  by 

Brms  ~  32  *^Torr  +  2  *GHz'  here  Erms  the  rms  electric  field 
in  volt/cm,  p  is  the  atmospheric  pressure  and  f  is  the  mi¬ 
crowave  frequency  (the  subscripts  indicate  the  units).  At  sea 
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level  (p  >>  f  ..  )  ,  the  microwave  power  flux  breakdown 
Torr  GHz 

2 

threshold  is  about  1,5  MW/ cm  .  This  kind  of  power  is  achiev¬ 
able  using  present  technology.  Free  electron  diffusion  is 
only  important  in  weakly  collisional  air  and  at  the  time  of 
breakdown,  and  can  be  ignored  once  the  plasma  builds  up;  this 
is  the  stage  we  are  most  interested  in. 

After  breakdown,  the  plasma  density  grows  exponentially 
in  time  quickly  until  the  density  is  high  enough  for  reflec¬ 
tion  and  absorption  of  the  microwaves.  The  microwave-plasma 
interaction  is  basically  dominated  by  electron-neutral  colli¬ 
sions  (to  be  described  in  Section  IV),  in  which  the  collision 

9 

frequency  is  roughly  given  by  vc  -  5  x  10  PTorr  at  suffi¬ 
ciently  high  microwave  power.  In  the  important  frequency 
range  of  a  few  GHz  to  several  hundred  GHz  the  system  is  very 
collisional  in  the  lower  atmosphere  (v  /to  >>  1,  where  to  =  2  irf 
is  the  angular  frequency  of  the  microwaves) .  In  this  case  the 
cut-off  plasma  density  of*  the  microwaves  is  given  by 

n/n_  -  v_/to,  here  n  *  1.26  x  10”8  f2  cm”3  is  the  critical 
c  c 

density,  instead  of  the  usual  n  ~  n  in  the  collisionless 

case. 

We  show  in  this  paper  that  reflection  depends  mainly  upon 
the  height  and  the  slope  of  the  density  profile.  In  the 
strongly  collisional  limit,  reflection  becomes  important  when 
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k_i  <  1  and  n/n  >  v  /w ,  here  k  is  the  wavenumber  of  the  mi- 

u  C  C  O 

crowaves  and  4  is  the  scale  length  of  the  density  gradient. 

In  the  weak,  collisional  limit,  they  are  given  by 

k  t  <  v/w  and  n/n  >  1.  In  useful  atmospheric  units,  the 
o  c  c 

conditions  for  the  two  limits  are  given  by  4(cm)  <  f ~^GHz  and 

"<”'3>  >  l0“  p,tmfGHz  vaUd  fot  patm  >  10~3  fGHz>  ana  t(cm)< 
10  Patm”1  ®nd  n(cm”3)  >  1010  f^jj2  valid  for 

Patm  <  10~3  fGHz  tespectively. 

In  vary  collisional  air,  the  condition  for  good  absorp¬ 
tion  is  given  by  (n/n_)  (k  A)  >  v _/w,  here  A  is  the  width  of 

c  c 

tbe  density  profile.  Rewritten  into  atmospheric  units,  this 
is  nA(cm”2)  >  1013  patm  or  ( n/N )  A  (cm)  >  10~6,  here  N  is  the 
atmospheric  density.  Therefore  only  low  ionization  of  the  air 
is  needed  to  give  good  absorption. 

The  currently  available  experimental  data  were  mainly  ob- 

* 

tained  in  dc  experiments  and  are  roughly  valid  in  the  range 

0  <  B  /P_  <  100  within  an  accuracy  of  10  to  20%. 

rms  Torr 

Checking  the  relation  of  the  data  to  the  kinetic  description 
carefully,  we  can  relate1  the  data  to  an  effective  microwave 
field  strength  which  is  defined  by  E(1  +  w2/v 

Therefore  the  rate  coefficients  (which  are  elastic  collision, 
ionization,  and  attachment)  and  thermal  energy  of  the  elec¬ 
trons,  can  be  determined  macroscopically  as  functions  of  mi- 
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crowave  field  strength  without  knowing  the  details  of  the  mi¬ 
croscopic  kinetic  description. 


We  present  hydrodynamic  calculations  for  both  planar  and 
spherical  geometries.  A  problem  that  must  be  avoided  is  the 
growth  of  the  plasma  density  in  all  space.  After  some  time 
the  density  even  grows  faster  at  places  close  to  the  microwave 
source  because  the  field  strength  is  stronger  there  due  to  re¬ 
flection.  Eventually  the  microwaves  will  be  cut-off  by  a  very 
high  plasma  density  near  the  source. 


In  the  planar  case,  we  present  two  types  of  calculations. 
First,  a  local  region  is  pre-ionized  by  other  means,  then  the 
pre-ionized  plasma  grows  to  high  enough  density  to  cause  high 
absorption  before  the  density  at  other  places  grows  signifi¬ 
cantly.  During  this  time,  a  steady  state  ionization  front  is 
formed  and  moves  rapidly  toward  the  microwave  source. 


*■ 


Second,  a  reflector  is  placed  at  a  boundary,  such  that 
standing  microwaves  are  formed  in  the  system.  The' incoming 
power  is  below  the  breakdown  threshold,  therefore  the  problem 
of  cutting-off  the  microwaves  near  the  source  can  be  avoided. 
However  the  maxima  of  the  standing  wave  can  be  above  the 
threshold,  and  the  multi-peak  plasma  density  is  formed  and  the 
microwaves  are  absorbed.  Stationary  state  absorption  is  ob- 
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tained  in  this  case. 

In  very  collisional  air,  results  for  both  cases  show  that 
the  product  n  (sum  for  multi-peak  density)  of  the  asymptotic 
plasma  density  profile  remains  roughly  the  same  at  fixed  mi¬ 
crowave  power  and  atmospheric  pressure,  and  is  quite  indepen¬ 
dent  of  other  parameters.  Therefore  it  gives  fixed  absorption 
(see  Eq.  49).  The  reason  is  that  the  microwaves  have  to  be 
attenuated  to  the  breakdown  value  Erms  ~  32  P<rorrin  the  densi¬ 
ty  profile  in  order  that  a  steady  state  is  achieved.  Thus  the 
required  attenuation  fixes  the  product  of  n  and  therefore  the 
absorption.  Energy  that  is  not  absorbed  will  be  transmitted 
in  the  pre-ionization  case  because  the  density  required  for 
significant  reflection  is  not  built  up. 

We  also  find  that  the  asymptotic  absorption  increases  ra¬ 
pidly  with  the  microwave  power  and  inversely  with  pressure. 
Total  absorption  occurs  for*power  not  very  much  higher  than 
the  threshold. 

In  weakly  collisional  air,  the  reflection  density  is 
easier  to  be  built  up.  We  also  find  that  the  reflection  de¬ 
creases  with  v^w,  the  collisionality  of  the  system,  and  is 
quite  independent  of  the  microwave  strength.  The  energy  that 
is  not  reflected  will  be  absorbed  due  to  the  longer  elongated 
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density  profile  at  the  back  in  the  pre-ionization  case. 

In  spherical  geometry,  which  simulates  converging  mi¬ 
crowaves,  the  cut-off  problem  is  also  avoided  because  the  mi¬ 
crowave  beams  are  .above  breakdown  threshold  only  near  the 
focal  region.  We  find  that  one  after  another  density  peaks 
are  created  at  the  end  close  to  the  focus,  and  they  move  to¬ 
ward  the  microwave  source.  However  the  density  peaks  cannot 
move  beyond  a  certain  point  in  the  system  because  the  field 
strength  is  below  the  threshold  beyond  it.  Asymptotic  absorp¬ 
tion  is  also  obtained  in  this  geometry,  however,  it  increases 
with  system  length  and  beam  convergence. 

In  all  planar  and  spherical  calculations,  the  electron 
temperature  tends  to  be  roughly  2  eV  in  steady  state  since  the 
microwave  strength  is  roughly  the  breakdown  strength. 

* 

* 

We  have  also  simulated  the  recent  focused  microwave  ex- 
2 

periments  by  two  connected  regions  with  different  spherical 
convergence  with  and  without  a  reflector  at  the  focus.  The 
ionization  occurs  initially  near  the  focus,  however  the  plasma 
grows  much  denser  near  the  threshold  region  at  later  time  and 
decouples  the  microwaves  from  the  focal  region.  The  length, 
location  and  motion  of  the  plasma  density  profiles  are  in  ex¬ 
cellent  agreement  with  those  of  the  experiments. 


The  organization  of  this  report  is  as  follows.  'The 
Boltzmann  equation  is  used  and  reduced  to  the  ‘pfoper  kinetic 
description  for  the  problem  in  Section  II.  In  Section  III  mo¬ 
ment  equations  are  obtained  for  the  problem  with  coefficients 
based  on  existing  experimental  data.  Arguments  are  made  to 
relate  the  data  to  the  microwave  field  strength.  In  Section 
IV  the  experimental  data  and  their  range  of  validity  are  re¬ 
viewed.  Reflection  and  absorption  of  microwaves  in  the  pres¬ 
ence  of  plasma,  treated  both  analytically  and  numerically,  are 
given  in  Sections  V  and  VI  respectively.  In  Sections  VII  and 
VIII,  microwave  propagation,  including  the  self-consistent  hy¬ 
drodynamics  of  the  plasma,  is  investigated  in  both  planar  and 
spherical  geometries.  Hydrodynamic  simulations  of  recent  fo¬ 
cused  microwave  experiments  are  described  in  Section  IX.  The 
work  is  summarized  in  Section  X. 
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XI.  DESCRIPTION  BY  THE  BOLTZMANN  EQUATION 

The  kinetic  theory  of  the  electrons  provides  the  essen¬ 
tial  description  of  the  microwave-atmospharic  interactions. 
The  high-frequency  phenomena  associated  with  the  motion  of  the 
Ions  and  neutrals  are  insignificant  due  to  much  higher  masses. 
Although  any  macroscopic  measurement  is  an  average  of  one  sort 
or  another,  the  microscopic  description  forms  the  basis  for 
the  understanding  of  the  measurements.  The  treatment  of  the 
electron  distribution  f  is  described  by  the  Boltzmann  equation 
which  is  given  by 

||  +  v  •  Vf  +  a  •  Vv  f  •  C  (1) 

where  v  and  a  are  the  velocity  and  acceleration.  The  colli- 
sional  term  C  represents  the  change  of  the  distribution  by 
elastic  and  inelastic  collisional  processes.  Elastic  colli¬ 
sions  mainly  result  in  the. transfer  of  momentum  and  inelastic 

V 

collisions  are  due  to  the  excitation,  ionization  and  attach¬ 
ment  of  the  electrons. 

The  process  of  simplifying  the  Boltzmann  equation  can  be 
based  upon  the  high  collisional ity  of  the  phenomena,  especial¬ 
ly  in  the  lower  atmosphere.  For  the  less  collisional  higher 
atmosphere,  the  limitations  of  the  description  should  be  care¬ 
fully  justified.  Since  collisions  take  place  in  times  of 
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order  of  10  second  (see  Section  IV)  in  the  lower  atmosphere, 

we  assume  that  any  substantial  ordering  of  the  direction  of 

the  motion  of  the  electrons  is  prevented  in  the  frequency 

range  considered,  and  the  velocity  distribution  is  nearly 

spherically  symmetrical  at  any  point  in  space.  Under  this 

consideration,  the*  electron  distribution  can  be  described  in 

3 

spherical  expansions 

£  »  fQ  +  ^  •  M  :  ?2  (2) 


where  f^  represents  the  spherically  symmetric  part,  and  the 
^2  and  ^2  represent  the  first  and  second  order  deviations  from 
spherical  symmetry,  and  <?  is  the  unit  velocity  vector. 
Similar  to  that  for  the  distribution  function,  the  expansion 
for  the  collision  term  is 
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by- 


*1 


*  «, 


JfT  *  v  +  **  H  X  ?!  +  I  V  v  .  f. 


(5) 


+  (v  *e  “  V  =  *i 


here  aB  is  the  acceleration  due  to  the  magnetic  field. 

Physically  Eq.  (4)  determines  the  balance  of  the  number 
density  and  energy  of  the  electrons.  Therefore  the  collision 
term  CQ  can  be  related  to  the  distribution  fQ  by  taking  ac¬ 
count  of  all  elastic  and  inelastic  collision  processes.  The 
elastic  collision  term  should  have  the  same  form  as  that  of 
the  third  term  in  equation  (4)  because  it  conserves  particles. 
Therefore  we  have3 

co  -  7  If  <v*  *<v>  fo>  -  V,  f„  (6) 

where:  A(v)  is  a  function  of  v  that  is  determined  fqr  elastic 

collisions  of  electrons  with  the  heavier  ions  and  neutrals. 

The  second  term  contains  v  which  relates  to  the  rates  of  ex- 

s 

citation,  ionization  and  attachment  of  electrons. 
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Eq.  (5)  determines  the  balance  of  the  electron  momentum. 
Therefore  the  collision  term  can  be  related  to  the 
electron-neutral  momentum  transfer  collision  frequency  vm  as 
follows 


Although  Eqs.  (4)  and  (5)  describe  only  the  zeroth  and 

the  first  order  distributions  of  the  spherical  expansion,  they 

are  exact  even  for  an  infinite  spherical  expansion3.  For  less 

collisional  phenomena,  in  which  the  ordering  expansion  fails, 

Eqs.  (4)  and  (5)  still  apply  if  the  phenomena  still  fits  the 

description  of  the  definition  of  the  distribution  given  by  Eq. 

-*• 

(2).  Right  now  there  is  no  ordering  among  tQf  and  ?2.  For 
both  collisional  and  collisionless  phenomena,  ?2  can  he 
modeled  physically  to  close  Eqs.  (4)  and  (5)* 

*v 

For  the  problem  of  interest,  we  assume  fQ  represents  only 

the  low  frequency  response  of  the  electrons.  This  is  the  same 

as  assuming  that  high  frequency  charge  separation  does  not 

♦ 

occur.  We  also  assume  ?2  is  only  the  low  frequency  response 
of  the  electrons  due  to  the  existence  of  the  high  frequency 
microwaves,  that  is  the  ponderomotive  effect  .  Finally  the 
high  frequency  response  due  to  the  microwaves  enters  only  in 


m 


V  V  ■*- 


•/  %•  <*  •: 


4 
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.  • 

** 

I,,  therefore  we  set 

*i  ■  V e  iut  +  V  <«> 

here  the  superscripts  h  and  L  represent  the  high  and  low  fre¬ 
quency  responses... 


Substituting  the  high  frequency  part  of  Eq.  (8)  into  Eq. 
(5)  and  assuming  there  is  no  dc  magnetic  field,  we  obtain  a 
linearized  high  frequency  equation  as  follows 


a_“  df* 

E _  _ o 

(id)  +  vj  dv 

in 


For  the  low  frequency  part,  Eq .  (5)  gives 


(9) 


+  v  Vf 

o 


■*  L 
+  «E 


ffo 

dv 


+ 


(10) 


We  have  neglected  the  second  to  last  term  in  Eq.  (5)  in  the 
above  equations  because  this  term  is  small  and  it  actually 
does  not  contribute  in  the  la^er  equations. 
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III.  MOMENT  EQUATIONS  FROM  THE  BOLTZMANN  EQUATION 


The  macroscopic  equations  can  be  deduced  from  the  kinetic 
equations  described  in  Section  II.  Quite  often  we  need  to 
refer  back  to  the  kinetic  theory  in  order  to  understand  cer¬ 
tain  macroscopic  measurements.  For  high  frequency  phenomena, 
we  also  need  an  equation,  obtained  from  Faraday's  and  Ampere's 
laws  of  Maxwell  equations,  to  describe  the  high  frequency 
electric  field  £h 


2„2  -+ 
C  V  Ei 


dire  ^  / 


v  f  dv 


(11) 


Here  we  have  set  V  •  Eh  s  0  since  high  frequency  charge 
separation  is  neglected  as  described  in  Section  II. 


Substituting  Eq.  (9)  into  Eq.  (11), 
microwave  propagation  is  obtained,  with 

kQ2  £  ?h  -  0 


the  equation  for 
*  eitot ,  as  follows 

(12) 


here  the  complex  dielectric  function  e  is  given  by 
e  -  1  -  Tr.-)dv 


dire2  .-Id  ,  v3 

3m  J  £o  "T  3v  in  -  iv 
v  m 


(13) 


In  general  v_  is  a  strong  function  of  v,  therefore  Eq.  (13) 
m 

is  difficult  to  evaluate.  However  we  will  use  Eq.  (13)  as  a 
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definition  for  an  effective  high  frequency  collision  coeffi¬ 
cient  vb  aa 

•  -  iv„  “JSft o  “I  Jv  (u  -  iv"_)d7  <14> 

n  v.  n 

here  n  s  /  f  dv.  Notice  that  when  the.  electron-neutral 
aoeentum  transfer  collision  frequency  vm  is  not  a  function  of 

v,  then  v.  *  v  . 

n  a 

Using  Eq.  (4),  the  particle  conservation  equation  can  be 
obtained,  after  integrating  over  all  velocity  space,  as  fol¬ 
lows 

|f  +  ^  •  ?  •  “  /  vs  f0  <*v  (15) 

here  J  *  1/3  /  v  dv.  The  source  term  is  due  to  ionization 
and  attachment  of  electrons. 

The  aomentum  conservation  equation  can  be  obtained  from 

*  * 

Eq.  (5)  *  after  multiplying"  by  mv/3  and  integrating  over  all 
velocity  space,  and  it  is  given  by 

m  +  Vp  +  en  EL  +  ^  •  P2  =  -  5  /  vmv  dv  ''  (16) 

2  ■*> 

here  a  is  the  electron  mass,  p  =  d/3  /  v  fQ  dv,  and 
P2  *2m/15  /  v  ?2  dv*  Here  p2  is  the  ponderomotive  pressure 

tensor,  which  cannot  be  obtained  in  our  present  derivation, 

4 

however,  it  can  be  derived  using  other  approaches  .  It  is 
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given  by* 

»  •  p2  -  (V  l^l2  -  2  Im  tuh  x  7  x  (vhS*)))  (17) 

where  a  1/3n  /  v  dv,  which  is  the  high  frequency  oscil¬ 
lation  velocity  ■  of  the  electrons  due  to  the  microwaves.  The 
second  term  on  the  right  hand  side  of  Eq.  (17)  can  be  larger 
than  the  first  term  if  vh  >>  “•  The  first  term  prevails  in 
the  presence  of  substantial  standing  microwaves. 

Basically  we  are  more  Interested  in  the  later  time  de¬ 
velopment  of  the  phenomena,  especially  as  the  system  ap¬ 
proaches  steady  state.  After  the  atmosphere  is  subjected  to 
microwave  breakdown,  electrons  and  ions  are  created  rapidly. 
As  soon  as  a  significant  number  of  ions  are  created,  the  mo¬ 
tion  of  the  electrons  is  limited  by  the  heavier  ions  through 
ambipolar  drag.  Therefore  we  set  j  =  nu,  here  u  is  the  hydro- 
dynamic  flow  velocity  of#fhe  ions.  Because  the  ions  do  not 
Interact  with  the  microwaves  the  hydrodynamic  description  for 
the  ions  is  appropriate.  The  momentum  equation  is  given  by 

M(-|^  nu  +  7  •  nuu)  -  en  a  0  (18) 

here  N  is  the  ion  mass,  and  the  ion-neutral  collision  and  the 
thermal  pressure  of  the  ions  have  been  neglected. 

*♦ 

Combining  Eqs.  (16)  and  (18),  E^  can  be  eliminated  and 
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m«  obtain 

•  ..  * 

M(|^  nu  +  7  •  nuu)  +7  •  P2  +  7P  +  5  ^  vm  v  ^  *  0 


The  evaluation  of  the  last  two  terms  cooes  from  the  ex¬ 
perimental  measurements  which  are  usually  done  in  the  steady 
state  and  in  absence  of  microwaves.  In  this  case,  a  moment 
equation  obtained  from  Eq.  (10)  is  given  by 


5  -  -  £  t. 


->  aS  -  |  v  /  fc  £  av 

i  m 


The  first  term  on  the  right  hand  side  of  Eq.  (20)  is  the 

mobility  term,  and  the  second  term  is  the  diffusion  term.  The 

c_d 

so-called  swarm  experiments  actually  determined  the  drift 
and  the  ratio  of  the  diffusion  coefficient  to  the  drift  inde¬ 
pendently  as  functions  of  an  applied  electric  field  tA> 
Therefore  an  effective  elect^gn  temperature  can  be  defined  as 

•A 

-  ■  '  fo  iZ/r  'oiR  <£>  <21> 

In  v  m 

as  a  function  of  Also  an  effective  dc  collision  frequency 
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The  best  way  to  approximate  the  last  two  terms  of  £q. 
(19)  is  to  relate  then  to  the  above  two  defined  quantities 
■easured  in  the. experiments,  and  Eq.  (19)  become 

M(j^  nu  +  V  •  nuu)  +  V  •  p2  +  9nee  +  ®vL  nu  =  0  (23) 

here  p^  is  given  by  Eq.  (17)  in  term  of  u^,  which  can  be 
derived  from  Eq.  (9)  to  be 


w  e 

Ti  m  (i<o  +  v^) 


(24) 


The  independent  variable  £  of  e  and  v_  can  be  set  to  be 

A  e  L 

equal  to  the  microwave  electric  field  as  (also  see  Section  I) 


/ 


1  +  “^h2 


(25) 


which  is  known  as  the  effective  electric  field.  The  meaning 
of  Eq.  (25)  is  that  £.  results  in  the  same  amount  of  colli, 
sional  energy  dissipation  in  both  dc  and  ac  cases,  therefore 
it  gives  roughly  the  same  ee  and  vL. 


Also  can  be  approximated  by  especially  in  the 
highly  eollisional  case.  Since  is  also  a  function  of  Eh, 
end  so  we  see  that  the  wave  equation  Eq.  (12)  is  a  nonlinear 
equation. 
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finally  the  continuity  equation  Eq.  (15)  can  be  written 

as 

t 

2  A 

+  V  •  nu  -  (vA  -  va)n 

here  v,  and  v,  are  the  ionization 
dents.  Both  werp  measured  in  dc 
of  the  dc  electric  field,  which  can  be  related  to  the  mi¬ 
crowave  field  defined  in  Eq.  (25)  using  the  above  argument. 


(26) 

and  attachement  rate  coeffi- 
9—13 

experiments  as  functions 
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IV.  RATE  COEFFICIENTS  FROM  EXPERIMENTS 

* 

in  order  to  investigate  the  dynamics  of  microwave  break¬ 
down  in  air#  experimentally  determined  rate  coefficients  are 

required.  Most  of  the  measurements  were  done  in  dc  expert- 

*  1 

■ants.  However#  there  are  ac  experiments  which  show  that  the 
breakdown  thresholds  of  air  is  similar  to  that  of  dc  experi¬ 
ments  at  high  pressures.  Within  the  accuracy  (10-20%)  of  the 
experiments#  several  approximate  formulae  are  modeled  from 
various  measurements  for.  the  rate  coefficients.  It  turns  out 
they  can  be  described  macroscopically#  without  knowing  the  de¬ 
tails  of  their  electron  velocity  dependence#  by  the  microwave 
field  strength. 

elastic  collision 

by  Eqs.  (14)  and  (22)  is 


(27) 


The 


el ect ron-ne  ut  r al 


frequency  vc 
given  by5'6 


*  vL  defined 


P. 


Torr 


5  x  109  /  ■ 

o  x  xu  .  0  +  g 


in  the  range  0  <  o  <  100.  Here  a  is  a  quantity  related  to  the 
effective  electric  field  described  in  Section  III  and  is  given 
by 


^Torr 


Brms _ 

/ 1  +«,2/vci 


(28) 


- - -  ,V, >>3 


J UXI 


1  •  n"  V  S./'  .* 

tlm  ”*n 


J Le-Xs 
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Motice  that  Bq .  (27)  is  a  nonlinear  equation  for  vc»  and  is 

roughly  constant  at  high  field. 

Once  the  electrons  are  heated  by  the  microwaves,  the 
processes  of  excitation,  ionization,  and  exchanging  energy 
with  the  ions  and  neutrals  occur.  Energy  balance  on  the  col- 
lisional  time  scale  will  determine  the  electron  temperature 
locally.  The  use  of  the  measured  temperature  includes  the 
electron  energy  gain  and  loss,  and  so  an  electron  energy  equa¬ 
tion  is  not  required.  The  electron  temperature  defined  by  Eq. 
(21)  is  given  by7,8 

ee3  -  2.1  x  l(f3  a  (91  +  a)  (29) 

in  the  range  0  <  a  <  40,  and  ee  is  in  eV.  The  temperature  is 
a  weak  function  of  a  as  expected  since  the  ionization  and  ex¬ 
citation  restrict  the  energetic  electrons.  Even  though  the 
range  is  not  large  we  an  extrapolate  to  larger  a  because  of 
the  weak  dependence  on  a.  *  * 

Ionization  occurs  only  when  the  microwave  breakdown 
strength  is  reached,  which  is  roughly'*’  a  «  32.  The  ionization 
rate  is  given  by3'1® 

^ -  -  5.14  x  1011  exp  (-73  a"0*44)  (30) 

pTorr 


in  the  range  32  <  o  <  100 
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Attachment  occurs  due  to  the  combination  of  electrons 
with  the  neutrals*  In  our  case,  tvre  processes  (e  +  02)  and  (e 
♦  202)  can  be  dominant.  The  two-body  attachment  rate  is  given 


11-13 


a2 

^Torr 


7.6  x  10"4  o2 (a  +  218) 2 


In  the  range  0  <  a  <  60  and  the  three-body  attachment  rate  is 

given  by13 


a3 

^Torr 


5  a 

5.7  x  10  p  ^m  273  173  (32) 

atm  a  +  6.1  a'*  +  1.6  a1'3  +  0.11 


in  the  range  0  <  a  <  5  which  corresponds  to  0  <  ee  <  1  eV. 
notice  that  is  only  important  at  very  low  microwave 
strengths.  Also  vft2  is  not  important  as  compared  to  the  ioni¬ 
sation  rate  for  larger  values  of  <*  after  the  breakdown. 


V.  REFLECTION  OF  MICROWAVES  BY  FIXED  ELECTRON 

t 

DENSITY  PROFILES 


The  reflection  of  microwaves  by  the  breakdown  plasma  in 
the  atmosphere  is  a  key  factor  in  the  determination  of  how 
much  microwave  power  can  be  delivered  through  the  atmosphere. 
Bare  we  will  calculate  the  reflection  of  the  microwaves  for 
various  electron  density  profles.  We  find  that  the  reflection 
basically  depends  on  the  height  and  the  slope  of  the  density. 
Both  analytical  and  manerical  calculations  will  be  presented. 


The  one-dimensional  wave  equation  (from  Eq.  12)  can  be 
written  as  (for  simplicity  the  subscript  is  dropped) 


d2E  .2 

+  ko  CE 


(33) 


°'  e  *  1  "  nc(l  -  ivcA>) 

Since  given  by  Eqs.  (27)  and  (28)  ,  is  a  function  of 

© 

Bf  therefore  the  wave  e qua t log  is  a  nonlinear  equation  of  E. 


When  the  reflection  is  very  weak,  v  can  be  considered  as 
a  constant  in  space  where  E  is  quite  constant  along  the  re¬ 
flection  path.  In  this  case  we  actually  can  analytically 
solve  Eq.  (33)  .  Two  types  of  plasma  density  profiles  are 
consiered.  The  first  type  is  an  exponential  density  profile 
which  is  given  by 


n 


exp  (x/f) 


(34) 
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The  second  typa  is  a  linear  density  profile  which  is  given  by 


Notice  that  l nidn/dxf1  *  L  n/n  is  the  local  density 

C 

dient  scale  length  for  both  profiles. 


(35) 

gra- 


The  reflection  can  be  obtained  analytically  for 

constant  for  the  above  density  profiles  by  using  the 

14 

Phase-Integral  Method  .  The  reflection  of  wave  energy  is 
given  by 

R  *  |  exp  (2ik0  dz)  T  (36) 


here  Z0  is  the  complex  point  where  the  complex  dielectric 
function  e  »  fl.  For  an  exponential  profile  Eq.  (35)  becomes1 

R  -  e’4  V  tan'1(v0/“>  ,37) 

14 

and  for  a  linear  profile  it  is  given  by 


R 


-  £  k  L  V  /u> 

»  3  o  \’C^ 


(38) 


From  Eqs.  (37)  and  (33),  it  can  concluded  that  reflec- 
tion  is  Important  when 

Jcoi  tan”1(vc/w)  <  1  (39) 

and 

n  v 

k  L  v  _/ui  «  k^£  ~  <  1 


(40) 
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In  fact,  Eq.  (40),  which  is  for  a  linear  profile,  is  roughly 
the  same  as  Eq.  (39)  if  vs  use  the  turning  point  density  (e  ■ 
0)  in  Eq.  (40).  From  Eq.  (33),  the  turning  point  density 
roughly  is 


(1  +  vj/o2)1'2 


(41) 


Using  the  turning  point  density  in  Eq..  (40),  we  find. 


k0i  v^w  (1  +  vcW)^2  <  1 


(42) 


One  can  see  that  Eqs.  (39)  and  (42)  are  roughly  comparable. 


In  general  Eq.  (39)  and  n  >  nfc  can  be  considered  to  be 

the  conditions  that  reflection  becomes  important.  In  the 

limit  v  /w  »  1,  they  are  given  by 
© 

kQi  <  1  and  n  >  ncvc/w  (43) 

and  in  the  limit  vc/w  <<  1,  they  are  given  by 

kQ»  <  ft)/vc  and  n  >  nc  (44) 

In  useful  atmospheric  units,  the  two  conditions  can  be  rewrit¬ 
ten  roughly  (valid  for  >  10~3  fGHz)  as 

£  (cm)  <  f'aGH3,  and  n(cm"3>  >  1013  PatmfG„2  <«) 

and  (valid  for  P  .  <  10“  3  fL„i  as 

atm  bnz 

f(cm)  <  l(f2  Patin_1  and  n(cm“3)  >  1010  f2GH2  (46) 
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Prom  Eqs.  (45)  and  (46),  one  can  see  that  reflection  is  im- 
portant  at  higher  atmospheric  altitude  and  lower  microwave 
frequency. 

The  above  analytical  calculations  are  only  valid  for  weak 
reflection.  For  a  complete  treatment  of  reflection,  numerical 
calculations  have  to  be  made.  Computations  are  made  for  a 
one-dimensional  finite  region.  The  boundary  conditions  are 
incoming  and  reflected  microwaves  at  the  left,  and  transmitted 
microwaves  at  the  right.  Eq.  (33)  is  solved  with  »C(E)  by 
iteration.  The  coefficients  of  reflection,  absorption  and 
.transmission  are  measured  at  the  boundaries. 

The  numerical  solution  in  Figure  1  shows  the  portion  of 

microwave  energy  reflected  just  before  reaching  a  point  on  the 

density  profile  as  a  function  of  the  local  density  at  that 

point  for  various  values  of^  and  L.  The  calculation  is  done 

in  very  collisional  air  (v_/u  ~  102)  with  p  *  1  atm,  f  •  3.2 

GHz,  and  ^mg/Bporr  “40.'  Obviously,  a  significant  portion  is 

reflected  after  n/n  <  v  /u  ,  which  is  the  second  condition  of 

c 

Eq.  (43)  . 

After  the  microwaves  have  passed  a  certain  density,  the 
reflection  reaches  its  maximum  value  (Figure  1)  which  is  a 
function  of  the  density  scale  length.  The  density  scale 


*.  *.  A  % 
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length'  for  the  maximum  value  is  shown  in  Figure  2.  The 

results  kQl  <  1  and.  kQL  <  10  2  (therefore  kQl  »  kQLn/nc<  1 

2 

since  n/nc  ~  ~  10  )  for  reflection  to  be  significant  are 

just  the  same  as  the  first  condition  of  Eq.  (43). 

Results  for  the  reflection  as  a  function  of  atmospheric 
pressure  (or  equivalently  as  a  function  of  microwave  frequen¬ 
cy)  are  presented  in  Figure  3.  Other  parameters  are  the  same 

2 

as  previously.  The  value  of  vc/w  is  roughly  10  Patm  in  this 

-2 

case,  therefore  the  region  <  10  corresponds  to  a  less 

collisional  atmosphere.  In  this  region  the  turning  point  is 
approximately  at  nc,  corresponding  to  the  second  condition  of 
Eq.  (44).  Also  L  *  1,  each  kQl  >  1  curves  for  both  profiles 
are  quite  similar  to  each  other  and  they  give  the  first  condi¬ 
tion  of  Eq.  (44)  . 

We  have  also  calculated  the  dependence  of  reflection  on 
'the  incoming  microwave  power.  It  was  found  that  the  depen¬ 
dence  is  very  weak  as  long  as  Ernis/pTorr  y>  8'  as  can  *een 
in  Eq.  (27).  This  is  the  usual  case  in  the  microwave  break¬ 
down  atmosphere. 
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VI.  ABSORPTION  OF  MICROWAVES  BY  FIXED  ELECTRON  DENSITY  PROFILES 

.  "  *  *  ‘  .  *  •  •  **  ' 

In  the  previous  section  we  have  calculated  the  reflection 
;of  microwaves  for  fixed  plasma  density  profiles.  Similar  cal- 
-  eolations  have  been  done  for  absorption.  Without  significant 
reflection,  the  absorption  can  be  represented  as  a 
phase- integral14  of  a  traveling  wave  which  is  given  by 


A  *  1  -  |  exp  (-i  kQ  /e  dx)  | 


In  general  Eq.  (47)  is  difficult  to  evaluate.  Here  we 

consider  a  simple  case  with  a  constant  plasma  density  slab  and 

constant  collision  frequency  v  .  The  absorption,  for  weak  re- 

c 

flection,  is  given  by  ^gj 

V  A  «  1  -  exp  k0A//(l  +  Vc2/w2)  (1  +  vc2/u>2  -  n/nj  ] 


y  C; 


here  A  is  the  width  of  the  sl»ab.  For  the  v  /w  >>  1,  the  ab- 

•  •  _  .  «  *!•  C 


sorption  is  given  by 


A  «  1  -  exp  (-  n/n_  w/v  k  A) 

c  CO 


The  highly  collisional  case  is  more  interesting  because 
it  gives  high  absorption.  For  comparison,  numerical  results 
(described  in  Section  V)  are  plotted  in  Figure  4  for  a  plasma 
slab  of  width  1  and  5  ,  here  X_  is  the  microwave  wavelength. 
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Other  parameters  are  p  s  1  atra,  f  =  3.2  GHz  and  6^/? yorr  s 

<10  (as  is  in  Figure  1  in  very  collisional  air)  .  One  can  see 

that  almost  total  absorption  is  obtained  before  the  turning 

2 

point  density  n/n_  10  ,  where  the  reflection  becomes  signi- 
ficant.  These  portions  of  the  absorption  curves  agree  very 
well  with  Eq.  (49)  if  we  use  an  average,  vc/u>  -  150.  At 
higher  density,  the  absorption  is  reduced  due  to  the  increase 
in  microwave  reflection.  However  the  reduction  is  less  if 
other  more  gentle  density  profiles  are  used,  since  here  the 
step  profile  gives  the  most  severe  reflection. 


If  the  reflection  is  not  serious,  Eq.  (49)  shows  that 
good  absorption  occurs  for 

5-V  <so> 

C 

Rewritten  in  atmospheric  units,  Eq.  (50)  becomes 

n  A  (cm’2)  >  1013  Patm  V*  (51) 

Since  the  atmospheric  density  is  N  3  x  1019  Pfttin  cm03,  there¬ 
fore  Eq.  (51)  becomes  "1^ 

g  A (cm)  >  10’6  (52) 

For  plasma  slabs  of  width  on  the  order  of  centimeters,  a  small 
percentage  of  ionization  gives  good  absorption. 


Although  only  a  constant  plasma  density  slab  is  used  here 
for  absorption  calculations,  other  density  profiles  will  give 
similar  results  with  an  effective  width  A  to  be  defined  as  the 
full  half-width  of  the  profile.  We  have  proved  this  numeri¬ 
cally. 


VII.  MICROWAVE  ABSORPTION  WITH  HYDRODYNAMICS  IN  PLANAR  GEOMETRY 

* 

All  the  previous  calculations  of  reflection,  absorption 
and  transmission  coefficients  use  only  fixed  plasma  density 
profiles.  In  a  realistic  situation,  the  density  profile 
should  be  determined  self-consistently  with  the  microwaves. 
The  time  evolution  of  the  density  profile  was  derived  in  Sec¬ 
tion  III,  and  the  one-dimensional  continuity  and  momentum 
aquations  are  given  by  Eqs.  (26)  and  (23)  as  follows 

+  ■§£  (nu)  *  “  va)  n  (53) 

and 

(54) 

M(|^  u  +  u  §£  u)  +  £  §£  n  ee  +  fp  +  mvcu  +  M<vi  "  va^u  **  0 
where  fp  is  the  ponderomot ive  term  which  is  given  by 

£p  *  J  Ki2  "  2  ST  ^“h  uh)J  (55) 

% 

* 

All  the  rate  coefficients  ,vc  and  ee  are  functions  of  E 

as  described  in  Section  IV.  Consequently  the  hydrodynamics  of 
the  system  evolves  according  to  the  microwave  fields  (Eq. 
33)  . 

Without  any  pre-arrangment ,  the  electron  density  usually 
grows  in  all  space  evenly  from  noise  once  the  microwave  power 
is  above  the  breakdown  threshold.  After  some  time,  the  densi- 


ty  tven  grows  faster  at  the  places  closer  to  the  microwave 
source  because  the  microwaves  are  stronger  there  due  to  re¬ 
flection.  Finally  the  microwaves  will  be  cut-off  due  to  the 
very  high  plasma  density  near  the  source.  This  is  not  the 
case  of  interest  because  after  a  while  the  microwaves  are 
blocked  from  leaving  the  source. 


the  important  issue  is  how  the  system  can  be  designed  to 
achieve  the  desired  physical  result.  Two  cases  will  be  inves¬ 
tigated  here.  First,  a  local  region  is  pre-ionized  by  other 

•  -  '  ,*■  « 

means,  and  second  a  reflector  is  placed  in  the  system.  The 

two  cases  will  be  examined  separately  in  the  following  subsec¬ 
tions. 


*3 
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K.  Pte-ionized  Localized  Plasma  Density 

If  we  Initialize  the  plasma  in  a  very  narrow  region  of 
relatively  high  density*  then  the  plasma  grows  to  high  enough 
density  to  cause  high  absorption  before  the  density  at  other 
places  grows  significantly.  Hopefully*  enough  energy  could  be 
absorbed  in  the  time  interval  before  the  plasma  grows  in  the 
regions  that  are  not  pre- ionized.  The  pre- ionization  could  be 
produced  by  other  devices. 

lie  will  start  with  strongly  collisional  cases.  The  first 

calculation  is  shown  in  Figure  5  with  ^zm^pTorr  *  40  (I  *  2.5 

2 

MW/ cm  )  *  p  »  1  atm*  f  -  3.2  GHz*  and  an  initially  localized 
plasma  density  n/nc  »  1  (n  -  1011  cm”3)  at  kQx  *  14  and  n/nc  - 
10  ^elsewhere.  The  position  of  the  pre-ionized  plasma  is  in¬ 
dicated  by  an  arrow  in  the  figure.  The  plasma  density  grows 
in  time  and  moves  toward  the  microwave  source.  Two  density 
profiles  .  at  later  times  t  -  50  nsec  and  100  nsec  are  shown  in 
the  figure.  This  motion  of  the  profile  is  mainly  due  to  the 
ionization  process  by  the  larger  microwave  field  in. front  and 
the  de- ionization  process  by  the  smaller  field  at  the  back  of 
the  profile. 

Roughly  starting  from  t  ■  100  nsec,  the  density  profile 
evolves  very  slowly  toward  steady  state  with  the  incoming  mi- 
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crowaves.  The  asymtotic  behavior  of  this  sort  can  be  seen  in 

the  tine  history  of  the  absorption  a3  shown  ih  Figure  6.  The 

asymptotic  maximum  density  is  ninax/nc*'  60  and  the  asymptotic 

absorption  is  about  0.82  with  little  reflection.  The  small 

reflection  is  consistent  with  the  prediction  in  Section  V 
2 

(here  v^w  -  10  )  .  At  t  -  200  nsec,  the  asymtotic  density  pro¬ 
file  is  destroyed  due  to  density  growth  evolved  at  locations 
other  than  the  preionized  region. 

The  notion  toward  the  source  is  very  rapid.  The  speed  is 
8 

about  10  cn/sec,  which  is  2  or  3  orders  higher  than  the  sound 
.  speed.  Therefore  it  is  due  to  the  rapid  ionization  rate.  The 
front  slope,  roughly  represented  by  an  exponential  profile 
with  k0&  -  1,  remains  about  the  same  for  all  time. 

.  Me  have  tried  to  find  how  the  initial  pre- ionization  af¬ 
fects  the  asymptotic  results.  It  turns  out  that  the  product 

.  .  *  *• 

of  the  naxinum  and  the  width  of  the  asymptotic  density  profile 
.-  and  the  absorption  remain,  about  the  same.  _  The  first  line  in 
Table  1  (a)  shows  the  results  for  the  case  we  just  described, 
and  the  results  for  an  identical  case  but  with  a  different  in¬ 
itial  pre-ionization  are  given  in  the  second  line.  The  same 
thing  happens  when  the  frequency  is  varied  as  shown  in  line  3 
in  Table  1  (a).  All  cases  give  the  roughly  same  nA  such  that 
the  value  of  jj—  ~  kQA  in  Eq.  (49)  is  found  to  be  roughly  the 
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same  and  therefore  the  absorption  is  the  same. 

The  above  constant  absorption  phenomena,  for  the  same  mi¬ 
crowave  power  and  atmospheric  pressure,  can  be  explained  as 
follows*  Near  the  peak  of  the  density  profile,  the  microwave 
field  strength  has  to  be  near  ~  32  PTorr  which  is  the 
breakdown  threshold,  so  that  steady  state  is  achieved.  The 
required  attenuation  of  the  incoming  microwave  field  to 
fixes  the  value  of  nA.  Therefore  fixed  absorption  is 
achieved  by  varying  either  the  maximum  or  the  width  of  the 
density  profile.  This  is  just  what  has  been  seen  in  the  hy¬ 
drodynamic  calculations. 


For  higher  microwave  power,  more  absorption  is  required 
to  attenuate  the  microwave  field  strength  to  E^.  Therefore 
higher  absorption  results  for  higher  power.  In  Table  1  (b) 

absorption  is  listed  for  various  microwave  powers.  One  should 
■’  note  that  absorption  increases  rapidly  until  the  microwaves 
are  totally  absorbed.  If  the  microwave  power  is  fixed,  ab¬ 
sorption  decreases  with  the  atmospheric  pressure.  The  depen¬ 
dency  goes  roughly  with  the  effective  microwave 


strength  o  «  E  /P„ 

*  rms  Torr 


We  also  find  that  the  density  required  for  significant 
reflection  (Eq.  (43))  is  not  attained  therefore  there  is 
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little  reflection,  and  that  the  density  scale  length  kQi  does 
not  vary  very  much  in  general  because  the  microwave  strength 

i  -f‘* 

is  attenuated  to  d  «  32  at  steady  state. 

Mow  we  turn  to  weakly  collisional  cases.  Since  the  con* 

•  • 

ditions  k  £  <  w/v_  and  n  >  n  (Eq.  44)  are  easier  to  meet, 
o  c  c 

here  reflection  dictates  the  microwave  absorption.  A3  before 
the  value  of  kQ£  does  not  vary  very  much  at  steady  state, 
therefore  the  reflection  depends  mainly  upon  v  /u  (Eq.  37)  or 

C 

the  collisional ity  of  the  system. 


The  parameters  used  in  the  hydrodynamic  calcalations  in 

the  weakly  collisional  cases  are  typically  I  *  33 

kW/cm2  (E  -  3500  volt/ cm) ,  f  »  35  GHz,  and  the  initial 

rms 

pre- ionized  density  is  0.1  nc  with  width  0.25  XQ.  In  this 

case  v  A)  »  17  P  .  .  The  calculations  show  that,  as  in  the 
c  atm 

collisional  cases,  the  plasma  density  front  moves  toward  the 
microwave  source  and  reaches  steady  state,  however  the,  profile 
at  the  back  continues  to  be  elongated  due  to  the  weak 
de-ionization,  rate  there.  Nevertheless  the  elongation  does 
not  affect  the  reflection  in  the  front  and  the  microwave  ener¬ 
gy  that  is  not  reflected  will  be  absorbed. 


Runs  with  various  pressures  (or  collisional ity)  are  shown 
in  Figure  7  in  which  the  weakly  collisional  regime  is  that  for 
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which  «.  0.1.  in  Figure  7  (a)  ,  the  asymptotic  density 

height  obtained  in  the  weakly  collisional  regime  is  roughly 

constant  since  a  «  0.783  E _ '  /f  „  does  not  depend  upon  pres- 

rms  ghz 

sure.  Also  the  asymtotic  absorption  increases  (reflection  de¬ 
creases)  with  pressure  as  predicted  in  the  previous  discus¬ 
sion.  The  cutoff  at  P  .  >  0.13  is  due  to  a  <  32  so  that  the 
poorer  is  below  the  breakdown  threshold.  The  speed  of  the  ion¬ 
ization  front  toward  the  microwave  source  is  shown  in  Figure  7 


v  v  / 
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B.  Reflector  at  the  Boundary 


If  a  reflector  is  put  at  the  right  boundary,  standing  mi¬ 
crowaves  will  be  formed.  The  incoming  power  is  below  the 
breakdown  threshold,  therefore  the  problem  of  cutting  off  the 
microwaves  near  the  source  mentioned  above  can  be  avoided. 
However  the  maxima  of  the  standing  wave  can  be  above  the 
breakdown  threshold,  so  that  plasma  is  generated  resulting  in 
absorption. 


A  eollisional  case  (Figure  3)  with  E  /P, 


rms  Torr 


20  (I 


'  0.62  MW/ cm  i  ,  P  -  1  atm  and  f  *  3.2  GHz  is  run  for  a  system  of 
length,  L  *  3.1  XQ.  Here  the  incoming  wave  a  *  20  is  below 
breakdown  threshold,  however  the  standing  wave  a  »  40  is  well 
'  above  threshold. 


Roughly  starting  from  tilths  t  »  100  nsec,  the  density  pro- 
file  evolves  very  slowly  *  towards  stationary  steady  state. 
Figure  8  shows  the  asymptotic  microwave  and  density  profiles 
at  t  *  200  nsec.  One  can  see  that  the  density  profile  forms 
with  spikes  aligned  with  the  anti-nodes  of  the  standing  wave. 
The  absorption  at  this  time  is  about  57%  such  that  the  maxima 
of  the  standing  wave  no  longer  is  o  *  49  but  o  -  32,  which  is 

the  breakdown  threshold,  so  that  a  steady  state  is  achieved. 


Calculations  also  have  been  performed  for  L  *  1.5  X  and 
6.2  *0.  The  asymtotic  results  for  the  three  cases  are  listed 
in  Table  2  (a)  .  •  One  can  see  that  the  product  of  the  length  of 
the  system  and  the  maximum  density  for  each  case  are  roughly 
the  same,  so  that  .the  absorption  is  roughly  the  same  (Eq. 
49).  This  is  the  same  constant  absorption  phenomenon  dis¬ 
cussed  in  the  previous  subsection.  In  this  case(  the  absorbed 
energy  is  distributed  among  the  density  spikes.  If  there  are 
a  large  number  of  density  spikes,  the  maximum  density  in  each 
spike  is  smaller. 

As  in  the  last  subsection,  the  absorption  increases  as 
the  power  increases.  This  is  listed  in  Table  2  (b)  ,  showing 
the  values  of  the  absorption  and  the  power.  Near  the  break¬ 
down  threshold,  the  microwaves  are  completely  absorbed. 

The  advantage  of  the  cast  with  reflector  is  that  there  is 

% 

absorption  of  the  microwaves  in  the  stationary  state.  Whereas 
in  the  previous  case  with  pre-ionization  the  moving  steady 
state  stays  only  for  a  finite  time. 


VIII.  MICROWAVE  ABSORPTION  WITH  HYDRODYNAMICS 
.  IN  SPHERICAL  GEOMETRY 

-The  problem  of  cutting-off  the  microwaves  near  the  source 

can  be  avoided  by  using  converging  microwave  beams  that  are 

above  breakdown  threshold  only  near  the  focal  region.  The 

-convergence  of  the  microwaves  can  be  obtained  from  structured 

or  several  antennas,  or  just  due  to  the  self-focusing  of  the 

microwaves.  The  actual  geometry  can  be  very  complicated,  for 

simplicity,  here  we  use  spherical  geometry  to  simulate  these 

effects.  We  assume  that  all  the  quantities  are  functions  of 

\t ,  the  spherical  coordinate,  only.  For  normally  incident  mi- 

15 

crowaves,  assuming  k  r  »  1#  Eq.  (12)  can  be  written  as 

o 

;  A-j  E  +  f  3?  E  +  kQ2eE  =  0  (56) 

Assuming  the  velocity  of  the  flow  is  in  r-direction  only, 
then  the  continuity  equatiorf  (Eq.  26)  becomes 

|£  +  (r2  nu)  =  <vA  -  va)n  (57) 

-  ■  T'  ^ 

The  momentum  equation  is  the  same  as  Eq.  (54)  in  the  previous 
section  with  x  replaced  by  r. 

We  have  run  a  case  with  Erms^Torr  3  30  (*  *  1*45 

2 

MW/ cm  which  is  just  below  breakdown  threshold) ,  p  *  1  atm  and 


M 
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£  *  3.2  GHz  in  a  system  that  extends  from  k  r  *  80  to  k  r  » 

100.  The  initial  plasma  density  is  n/n,  *  1  everywhere  (this 

c  % 

does  not  affect  the  later  stage)  .  The  density  and  microwave 

profiles  are  shown  in  Figure  3  at  time  t  ■  100  and  500  nsec 

~  respectively.  The  density  grows  initially  at  the  far  right 

until  transmission  of  the  microwaves  is  blocked  to  a  value 

Just  below  the  breakdown  threshold.  The  density  peak  grows  to 

-  V  n/n„  ~  7  and  moves  toward  the  microwave  source.  At  t  »  500 
.  c 

/'./  nsec  the  density  peak  moves  very  close  and  slowly  to  the 
source  such  that  the  transmitted  microwaves  pick  up  some 
:  strength  at  the  far  right  again  to  initialize  another  density 
peak.  The  process  repeats  itself.  At  this  time  the  absorp- 
tion  oscillates  about  an  asymptotic  value  of  about  18%  with 
little  reflection. 


Another  system  with  60  <.  kQr  <  100,  which  is  twice  as 


long,  has  been  run  with  the  same  parameters.  Figure  10  shows 


the  late  time  profiles  at  \  »  1250  nsec.  Many  density  peaks 
have  emerged  and  move  toward  the  source.  The  density  peak 
cannot  move  beyond  kQr  9  90  because  the  microwave  strength  is 
below  the  breakdown  threshold  to  the  right  of  that  point.  The 
microwaves  in  the  system  are  attenuated  to  the  breakdown 
value.  The  absorption  oscillates  about  an  asymptotic  value  of 
about  56%  with  little  reflection. 


In  Table  3  the  asymptotic  absorption  in  systems  of  dif¬ 
ferent  lengths  is  listed.  We  see  that  the  longer  the  system 
the  higher  the  absorption,  and  absorption  increases  as  the  mi¬ 
crowave  convergence  increases  (line  4). 
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IX.  SIMULATIONS  OP  FOCUSED  MICROWAVE  EXPERIMENTS 

> 

In  the  experiments  'carried  out  at  NRL  (Figure  11),  the 
microwaves  (frequency  35  GHz  and  Power  112  kW)  are  focusd 
through  a  lens  (7*6  cm  diameter  and  11.2  cm  focal  length)  into 
the  nitrogen  gas  (25  Torr) •  The  measured  elliptic  area  of  the 
focus  is  v  x  0.55  cm  x  0.75  cm.  Using  a  circular  area  to 
model  the  measured  elliptic  area,  the  equivalent  radius  is 
0.64  cm.  The  radial  profile  of  the  microwaves  is  approximated 
by  a  Gaussian  with  an  e-folding  distance  of  0.77  cm.  The  max¬ 
imum  power  flux  is  60  kW/cm  (equivalent  to  the  total  power 
evenly  distributed  in  a  circle  of  radius  0.77  cm)  at  the  axis. 

The  diffraction  of  the  microwaves  starts  to  occur  at  an 
axial  distance  of  3.5  cm  (4  wavelengths)  from  the  focus. 
Between  this  point  and  the  lens,  the  microwaves  are  focused 
spherically.  Our  simulation  starts  at  an  axial  distance  of 

•bout  6  cm  (7  wavelengths)  from  the  focus  with  I  »  8.2 

2 

kW/cm  at  the  axis  (equivalent  to  the  total  power  evenly  dis¬ 
tributed  in  a  circle  of  radius  2.1  cm).  The  microwaves  are 
focused  spherically  in  a  distance  of  2.5  cm  and  connected  to 
the  solution  in  the  diffraction  region  (3.5  -*■  0  cm)  .  We  ap¬ 
proximate  the  diffraction  effect  by  another  spherical  effect 
of  aspect  ratio  1.55  (equivalent  to  a  9.77  cm  radius  circle  at 
focus  and  1.2  cm  radius  circle  at  3.5  cm)  .  The  schematic  of 
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the  simulation  region  is  shown  in  Figure  11. 

Although  the  coefficients  in  our  computer  code  are  for 
air,  the  code  can  be  used  to  approximate  N  gas.  As  in  the 
experiments,  two. cases  are  run  with  and  without  a  metalic  re¬ 
flector  placed  at  the  focus.  The  microwave  pulse  length  is 
,  1  Msec. 


In  the  case  without  reflector,  the  spatial  profiles  are 
presented  in  Figures  12  and  13  (microwaves  incoming  from  the 
left)  •  At  t  <■  50  nsec,  the  ionization  of  the  air  by  the  mi¬ 
crowaves  occurs  mainly  near  the. focus  (Pigure  12  a),  however, 
at  the  same  time  the  microwaves  are  attenuated  there  due  to 


the  presence  of  the  plasma  (Figure  13  a)  •  The  associated 
electron  temperature  is  plotted  in  Figure  13  (b)  • 


At  t  »  100  nsec,  the  p^fesma  grows  faster  near  the  connec- 

■*! 

tion  region  (Figure  12  b)  because  the  microwaves  are  stronger 
there.  At  t  ■  400  nsec,  the  plasma  at  the  ionization  front 
reaches  the  collisionless  critical  density  and  becomes  thicker 
(Figure  12  cl  ,  thus  the  microwaves  are  decoupled  from  the 
focal  region.  The  ionization  front  moves  slowly  toward  the 
microwave  source  until  at  t  ■  1  ysec  (Figure  12  di  ,  and  the 
front  cannot  move  beyond  6  cm  because  the  microwave  strength 
is  below  the  breakdown  threshold  beyond  that  point.  The  de- 
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coupling  of  the  microwaves  from  the  focal  region  at  1  psec  is 
shown  in  Figure  13  (a)  . 

The  calculated  profiles  agree  quite  well  with  those  from 
the  experiments  without  a  reflector  (the  frame  pictures 
presented  in  Figure  2  of  reference  2)  .  The  common  observa¬ 
tions  are  (1)  the  initial  less  dense  plasma  near  the  focal  re¬ 
gion,  (2)  the  much  denser  plasma  just  inside  the  breakdown  re¬ 
gion  at  later  times,  (3)  the  slow  motion  of  the  denser  plasma 
at  later  times,  and  (4)  the  length  of  the  ionization  region 
being  about  5  cm. 


Our  calculated  plasma  density  (4  x  IB12  -  2  x  IB13  cm**3) 
is  also  within  the  experimental  measured  range.  However,  our 
calculated  electron  temperature  is  roughly  an  eV  lower.  The 
higher  experimental  value  is  probably  due  to  the  use  of  15 
Torr  Nj  and  19  Torr  He  at  the  time  of  the  measurement  of  the 
electron  temperature.  ** 


The  calculated  spatial  profiles  are  presented  in  Figures 
14  and  15  for  the  case  of  a  metal  reflector  at  the  focus  per¬ 
pendicular  to  the  incoming  microwaves.  Because  of  the  stand¬ 
ing  waves  set  up  by  the  reflector,  the  initial  ionization  is 
stronger  and  spatially  more  non-uniform  than  in  the  previous 
case.  This  is  shown  in  Figure  14  (a)  for  an  early  time,  t  ■ 


10  nsec.  Die  standing  wave  characteristic  is  shown  in  Figure 

15  (a)  with  the  microwaves  attenuated  near  the  £ocus  due  to 

1  >•*' 

the  presence  of  the  plasma. 

As  in  previous  cases,  the  plasma  grows  faster  near  the 
connection  regioh  and  moves  out  slowly  as  shown  in  Figure  14 
(b)  ,  (c)  ,  and  (d)  at  t  -  100,  400  and  1000  nsec  respectively. 
However,  the  denser  plasma  profile  in  this  case  is  split  due 
to  the  growth  of  the  spikes.  The  decoupling  of  the  microwaves 
from  the  focal  region  at  t  *  1  ysec  is  shown  in  Figure  15  (bl  . 

...»  ?  ■  •  .. 

Qualitatively  the  cases  without  and  with  reflector  are 
similar.  This  is  consistent  with  the  observations  in  the  ex¬ 
periment  (Figures  2  and  3  in  Reference  2  for  the  two  cases)  • 
The  additional  features  due  to  the  presence  of  the  reflector 
.shown  in  the  frame  pictures  are  (1)  stronger  ionization  in  the 

focal  region  at  early  time  and  (2)  the  spiky  feature.  In  gen- 

V 

oral,  similar  results  are  obtained  at  later  times  due  to  the 
strong  absorption  (A  90%)  of  the  microwaves  by  the  plasma, 
such  that  the  reflection  by  the  reflector  is  not  important. 

Die  motion  of  the  ionization  front  and  the  decoupling  of 
the  microwaves  can  be  seen  in  the  time  plot  of  the  position  of 
the  maximum  of  the  electric  field  in  Figure  16  (similar 
feature  for  the  case  without  reflector)  .  As  quickly  as  59 


nsec,  the  microwaves  are  decoupled  from  the  focal  region. 
After  that  the  ionization  front  moves  slowly  just  inside  the 
breakdown  region.  ‘Hie  speed  of  the  ionization  iront  is  obvi¬ 
ously  a  function  of  time.  The  average  speed  is  4  to  5  x 
10**  cm/sec.  The  time  plots  agrees  with  the  experimental  ob¬ 
servations  (Fig ore  3  b  in  Reference  2). 
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SUMMARY 


In  this  report  the  air  breakdown  due  to  microwaves  propa¬ 
gating  in  the  atmosphere  is  investigated.  The  kinetic  and  hy¬ 
drodynamic  descriptions  of  the  problem  are  reduced  from  the 
Boltzmann  equation.  The  available  experimental  data  for  rate 
coefficients  are  used  and  related  to  the  microwave  field 
strength  macroscopically.  The  important  processes  are  elec¬ 
tron-neutral  collisions,  ionization,  attachment,  and  excita¬ 
tion..  Thermal  equilibrium  of  the  electrons  with  the  mi¬ 
crowaves  is  also  obtained. 

Reflection  and  absorption  of  the  microwaves  in  the  pres¬ 
ence  of  the  plasma  density  profiles  are  investigated  both  an¬ 
alytically  and  numerically.  Conditions  for  reflection  and  ab¬ 
sorption  are  obtained. 

. 

Hydrodynamic  calculations  in  planar  geometry  were  carried 
out  for  two  cases:  pre-ionization  in  a  locallized  region  ini¬ 
tially  and  with  a  reflector  at  a  boundary.  In  the  first  case, 
enough  microwave  energy  can  be  absorbed  in  the  time  interval 
of  interest,  whereas  a  stationary  state  can  be  achieved  in  the 
second  case.  Calculations  in  spherical  geometry  show  that 
■any  density  peaks  emerge  from  the  end  close  to  the  focus  and 
tfeccelerate  toward  the  microwave  source. 


Some  typical  planar  results  show  that  the  air  is  ionized 
rapidly  to  form  an  ionization  front  moving  toward  tha  mi¬ 
crowave  source  which  consequently  decouples  the  microwaves 
from  the  original  region.  By  focusing  the  microwaves  or  plac¬ 
ing  a  reflector,  localized  regions  can  be  heated  by  the  mi¬ 
crowaves. 

The  electron  temperature  obtained  tends  to  be  roughly  2 

•V  at  steady  state  regardless  of  the  microwave  power.  The 

heated  electrons  are  able,  to  exchange  energy  with  the  neutral 

16 

atmospheric  molecules  and  will  possibly  enhance  the  neutral 
temperature. 

At  high  altitude,  the  heating  of  electrons  by  other 

mechanisms,  such  as  the  collective  plasma  effects,  is  possi- 

17 

ble.  TWo  effects,  parametric  instabilities  and  resonant  ab- 
18 

sorption  ,  can  result  in  heating  suprathermal  electrons  with 

‘r*  + 

higher  temperature.  These  effects  are  being  investigated. 

The  limitations  of  the  calculations  are  due  to  the  avai¬ 
lability  of  the  experimental  data  for  the  rate  coeffients. 
Nevertheless,  comparison  of  the  calculations  with  recent  mi¬ 
crowave  experiments  give  excellent  agreement.  More  accurate 

19 

rate  coefficients  are  being  Investigated  in  the  presence  of 
a  microwave  discharge. 
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Table  1 


(a) 

f  (GHz) 

n/nc 

A/Xe 

absorption 

3.2 

60 

0.6 

0.82 

3.2 

27 

1.3 

0.80 

32.0 

t  0.6 

0.82 

E  /P 

rms'  Torr 

...  Power  (MW/cm2 ) 

absorption 

35 

1.9 

0.33 

40 

’  ■  2.5 ;  -• 

0.82 

45 

3.1 

0.97 

50 

■  ;  3.8 

0.99 

60 

*'5.5 

1.00 

Asymptotic  values  of  density  height,  width  and  microwave  ab¬ 
sorption  from  planar  hydrodynamic  calculations  in  very  colli- 
aional  air  with  pre-ionization,  (a)  Results  obtained,  lines  1 
and  2,  by  varying  initial  pre-ionization  and,  lines  1  and  3, 
by  varying  frequency.  They  all  give  roughly  ths  same  nA  and 
absorption.  Here  p  a  1  atm  and  I  =  2.5  MW/ cm2,  (b)  Absorp¬ 
tion  increases  as  power  increases.  Here  p  =  1  atm  and  f  =  3.2 


Table  2 


4? 

Page ,52^ 


absorption 


0.60 


0.57 

0.62 


S  /P 
nas'  Torr 


Power  (MW/cm  )  absorption 


0.62 


0.60 


0.96 

1.40 


0.89 


0.99 


Asymptotic  values  of  density  height  and  oicrowave  absorption 
fro*  planar  hydrodynamic  calculations  in  very  collisional  air 
with  reflector,  (a)  Results  obtained  by  varying  the  system 
length.  They  all  give  roughly  the  same  nL  and  absorption. 
Here  p  s  1  atm,  f  =  3.2  GHz  and  I  =  0.62  Mtf/cm2.  (b)  Absorp¬ 
tion  increases  as  power  Increases.  Here  p  and  f  are  the  same 

% 

as  in  (a)  . 
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Table  3 


Region  o£  computation 


80  <  k  r  <  100 

—  o  — 

60  <  k  r  <  100 

—  o  — 

40  <  k  r  <  100 

—  o  — 

30  <  k  r  <  50 

—  o  “ 


absorption 


0.18 


0.56 

0.80 


0.55 


The  asymptotic  values  of  microwave  absorption  from  -spherical 
hydrodynamic  calculations  for  systems  of  different  lengths. 
The  longer  or  the  more  convergent  the  system  is  the  higher  the 
absorption.  Here  p  ~  1  atm,  f  =  3»2  GHz  and  I  s  1.*I5  MW/ cm  . 
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The  portion  of  microwave  energy  reflected  just  before  reaching 
a  point  on  the  plasma  density  profile  as  a  function  of  the 
local  density  at  that  point.  The  curves  are  for  exponential 
(with  1),  linear  (with  L) ,  or  step  profiles.  The  calculation 
Is  done  in  very  oollisional  air  with  p  =  1  atm,  f  s  3*2  GHz, 
and  I  .  2.5  MW/ca2  (E^g/p^^  =  40). 


S  .V  .  . 


1*1 


?  I0"4 


I»- 


Page 


0.4  0.6 

Reflection 


J|0-3 

i.O 


FIGURE  2 


The  maximum  reflection  from  Figure  1  as  a  function  of  plasma 
density  scale  length. 
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FIGURE  3 

•  •  • 

Reflection  as  a  function  of  pressure  for  various  plasma  densi¬ 
ty  scale  lengths.  This  case  has  the  same  f  and  E/P  as  in  Fig¬ 
ure  1.  The  boundary  between  strong  and  weak  collisionality  is 

_2 

roughly  at  p  <  10  ata. 


FIGURE  4 

Absorption,  transmission  and  reflection  of  microwaves  as  func¬ 
tions  of  maximum  density  for  a  plasma  slab  of  width  and 
5  XQ.  This  is  a  very  collisional  case  that  has  the  same  p,  f 
and  I  as  in  Figure  1. 


FIGURE  5 


Density  and  alerowave  profiles  at  different  times  (in  nsec). 
The  initial  pre-ionized  localized  density  n/nc  s  1  is  indicat¬ 
ed  by  an  arrow,  and  it  evolves  to  an  asymptotic  density  pro- 
# 

file  (t  =  100)  that  moves  rapidly  toward  the  microwave  source. 
This  is  a  very  collisional  case  that  has  the  same  p,  f  and  I 
as  in  Figure  1* 
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FIGURE  7 

Asyaptotic  values  obtained  in  runs  with  f  =  35  GHz,  I  s  33 

KW/ca2,  and  initially  pre-ionized  density  n  a  0.1  n  with 

c 

width  0.25  \Q.  The  systea  is  weakly  collisional  for 
pat»  <<  ®#1*  Displayed  are:  (a)  Maxmimum  density  and  absorp¬ 
tion,  and  (b)  speed  of  the  ionization  front  as  functions  of 


pressure 
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FIGURE  8 


Density  end  microwave  profiles  at  t«200  nsec.  A  stationary 
state  is  achieved  by  using  a  reflector  at  the  right  such  that 
only  the  initial  maxima  of  the  standing  microwaves  are  above 
the  breakdown  threshold.  This  is  a  very  collisional  case  with 
pal  atm,  f  s  3.2  GHz  and  I  s  0.62  MW/ cm2. 
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FIGURE  9 


Density  and  microwave  profiles  at  two  different  times  (in 
nsec)  for  spherical  geometry.  Density  peak  is  created  at  far 
right  and  moves  to  near  the  left  edge  where  the  microwave  flux 
la  below  the  breakdown  threshold.  This  is  a  very  collisional 
ease  with  p  s  1  atm,  f  s  3.2  GHz  and  I  s  1.45  AM cm2. 
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FIGURE  10 


Density  and  microwave  profiles  at  t  :  1250  nsec  for  the  system 
twice  as  long  as  that  in  Fig.  9*  Many  density  peaks  are 
created  at  far  right  and  bunch  up  near  the  left  edge.  Note 
that  the  microwaves  are  attenuated  to  the  breakdown  value. 
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FIGURE  11 

Schematic  of  the  geometry  for  calculations  to  simulate  the  fo¬ 
cused  microwave  experiments  with  p  s  25  Torr,  f  :  35  GHz  and 
power  112  kW.  The  estimated  power  fluxes  on  the  axis  are 
shown  and  diffraction  (simulated  by  another  spherical  effect) 
starts  roughly  3*5  cm  from  the  focus. 


FIGURE  12  • 


Density  profiles  at  $  different  times  for  the  case  without  re¬ 
flector.  The  ionization  occurs  initially  near  the  focus, 
however,  the  focal  region  is  decoupled  from  the  microwaves  due 
to  the  stronger  ionization  near  the  connection  region  at  later 
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*'  FIGURE  13 

Microwave  and  temperature  profiles  corresponding  to  the  plots 
shown  In  Figure  12.  Microwaves  are  decoupled  from  the  focal 
region  as  shown  at  t  a  1000  nsec. 
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Density  profiles  at  4  different  times  for  the  case  with  a  re 
flector  at  focus.  A  spiky  structure  is  observed  in  this  case 
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FIGURE  16 

Time  history  of  the  position  of  the  maximum  of  the  microwave 
electric  field.  This  shows  the  microwaves  are  decoupled 


quickly  at  about  50  nsec  and  the  ionization  front  moves  very 
slowly  near  the  connection  region. 
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